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This  paper  presents  a  review  on  the  state-of-the-art  maximum  power  point  tracking  (MPPT)  techniques 
for  PV  power  system  applications.  The  main  techniques  that  will  be  deliberated  are  the  Perturb  and 
Observe,  Incremental  Conductance  and  Hill  Climbing.  The  coverage  will  also  encompass  their  variations 
and  adaptive  forms.  In  addition,  the  more  recent  MPPT  approaches  using  soft  computing  methods  such 
as  Fuzzy  Logic  Control,  Artificial  Neural  Network  and  Evolutionary  Algorithms  are  included.  Whilst  the 
paper  provides  as  thorough  treatment  of  MPPT  at  normal  (uniform)  insolation,  its  focus  will  be  on  the 
applications  of  the  abovementioned  techniques  during  partial  shading  conditions.  It  is  envisaged  that 
this  review  work  will  be  a  source  of  valuable  information  for  PV  professionals  to  keep  abreast  with  the 
latest  progress  in  this  area,  as  well  as  for  new  researchers  to  get  started  on  MPPT. 
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1.  Introduction 

Despite  the  relatively  high  cost  of  solar  modules,  PV  power 
systems  have  been  commercialized  in  many  countries  for  its  long 
term  economic  prospects  and  more  crucially,  the  concerns  over 
the  environment  [1],  These  systems,  which  could  be  grid- 
connected  or  stand-alone,  are  being  installed  in  wide-ranging 
power  capacities  and  by  using  various  silicon  technologies. 
Regardless  of  the  type  and  size,  one  critical  component  of  any 
PV  is  the  effectiveness  of  its  maximum  power  point  tracker 
(MPPT).  This  area  has  been  and  is  still  attracting  immense  interest 
from  PV  research  communities  as  well  as  industrial  players 
because  it  is  the  most  economical  way  to  improve  the  overall 
PV  system  efficiency. 

Conceptually,  MPPT  is  a  simple  problem — it  is  basically  an 
operating  point  matching  between  the  PV  array  and  power 
converter.  However,  because  of  the  non-linear  I-V  characteristics 
of  the  PV  curve  and  the  consequence  of  the  varying  environ¬ 
mental  conditions  (particularly  insolation  and  temperature), 
tracking  the  correct  maximum  power  point  (MPP)  can  sometimes 
be  a  challenging  task.  The  tracking  becomes  more  complicated 
when  the  entire  PV  array  does  not  receive  uniform  insolation — a 
condition  known  as  partial  shading  [2].  This  issue  has  received 
much  attention  lately,  partly  due  to  the  proliferations  of  building- 
integrated  (B1PV)  in  urban  areas.  The  B1PV,  which  is  commonly 
installed  on  rooftop,  often  lends  itself  to  shading  because  of 
unfavourable  physical  constraints  such  as  space  limitations  and 
unforeseen  structural  development.  Notwithstanding,  partial 
shading  is  also  being  experienced  in  large  PVPS  installations, 
i.e.,  solar  farms.  Typically,  it  is  caused  by  the  clouds  that  strike  on 
certain  spots  of  the  solar  array,  while  other  parts  are  left 
uniformly  irradiated.  Such  meteorological  condition  is  quite 
common  in  tropical  and  near-equatorial  regions.  Another  source 
of  partial  shading-like  characteristics  is  exhibited  by  module 
irregularities;  a  common  example  would  be  the  presence  of 
cracks  on  one  or  more  modules  of  the  PV  array. 

When  a  PV  system  is  subjected  to  partial  shading,  its  P-V 
curves  exhibit  multiple  peaks  with  several  local  peaks  and  a 
unique  global  peak  (GP).  The  conventional  MPPT  algorithm,  such 
as  Perturb  and  Observe,  Hill  Climbing  and  Incremental  Conduc¬ 
tance  (to  name  a  few)  is  likely  to  be  trapped  at  one  of  the  local 
peak,  simply  because  its  algorithm  could  not  differentiate  the 
local  with  the  GP.  Consequently,  it  oscillates  around  the  vicinity  of 
the  local  peak  and  remains  there  indefinitely — a  significant 
reduction  in  the  PV  power  yield  will  be  experienced  [3]. 

Over  the  years,  numerous  MPPT  techniques  have  been  pro¬ 
posed;  this  is  evident  from  the  growth  in  number  of  related  papers 
published  in  various  journals  and  proceedings  [4,5].  Despite  the 
fact  that  these  methods  are  designed  for  the  same  objectives,  they 
differ  markedly  in  terms  of  complexity,  convergence  speed,  steady 
state  oscillations,  cost,  effectiveness  and  flexibility.  Furthermore, 
numerous  researches  have  been  carried  out  to  address  the  partial 


shading  effect  by  improvising  the  conventional  MPPT  controllers. 
However,  with  recent  availability  of  vast  and  low  cost  computing 
power,  new  methods  based  on  soft  computing  techniques  have 
been  attracting  considerable  interests.  This  approach,  which 
includes  the  Fuzzy  Logic  Controller,  Artificial  Neural  Network  and 
Evolutional  Algorithm,  offers  numerous  opportunities  for  more 
robust  and  flexible  MPPT  schemes,  in  a  way  that  is  not  possible 
using  conventional  techniques.  More  importantly,  these  methods 
are  well  suited  to  cater  for  the  partial  shading  conditions. 

In  a  short  span  of  time,  the  body  of  knowledge  related  to  this 
subject  has  grown  tremendously.  Clearly,  it  is  difficult  to  keep 
track  with  the  literatures  unless  a  single  reference  with  concise 
and  accurate  summaries  of  the  various  MPPT  schemes  is  made 
available.  The  most  recent  reviews  [4,5]  are  carried  out  over  five 
years  ago,  and  obviously  progress  since  then  has  been  consider¬ 
able.  Therefore,  it  seems  appropriate  at  this  juncture  to  provide  a 
comprehensive  and  up  to  date  review  of  the  subject.  In  particular, 
the  MPPT  for  partial  shading  conditions  is  not  fully  documented; 
justifiably,  there  is  a  critical  need  to  carry  out  this  exercise 
even  more. 

Although  this  paper  will  mainly  deal  with  the  latest  work  on 
MPPT,  it  is  envisaged  that  some  background  knowledge  would  be 
helpful  for  certain  group  of  readers  to  follow  through  this  paper 
effectively.  Accordingly,  information  on  the  principles  of  opera¬ 
tion,  limitations  and  advantages  of  the  conventional  MPPT  meth¬ 
ods  (and  their  variations)  are  provided.  However,  since  the 
number  of  papers  related  to  one  particular  method  is  quite  large, 
it  is  imperative  that  only  works  with  significant  contributions  are 
cited.  Papers  that  refer  to  previous  works  with  minor  modifica¬ 
tions  or  improvements  may  not  be  included  in  the  reference  list. 
In  that  regard,  apologies  are  offered  to  the  respective  authors. 


2.  MPPT  control  structure  with  power  converters 

The  aim  of  employing  MPPT  is  to  ensure  that  at  any  environ¬ 
mental  condition  (particularly  solar  insolation  and  temperature), 
maximum  power  is  extracted  from  the  PV  modules.  This  is 
achieved  by  matching  the  PV’s  MPP  with  the  corresponding 
power  converter’s  operating  voltage  and  current.  Fig.  1  shows 
the  general  block  diagram  of  MPPT  in  conjunction  with  a  dc-dc 
converter.  Although  a  stand-alone  dc-dc  system  is  depicted  here, 
the  application  can  be  extended  to  a  grid  connected  PV  system  by 
adding  other  power  electronic  devices  such  as  inverter  and  grid 
components.  Basically,  the  MPPT  works  by  sensing  the  current 
and  voltage  of  the  PV  array;  using  this  information  the  array 
power  is  calculated  and  compared  with  the  present  value  of  MPP. 
Accordingly,  the  duty  cycle  of  the  converter  is  adjusted  using  a  PI 
or  hysteresis  controller  to  match  the  MPP — which  in  turn,  forces 
the  converter  to  extract  the  maximum  power  from  the  array. 

An  alternative  control  structure  is  characterized  by  directly 
updating  the  duty  cycle  of  the  power  converter;  this  is  known  in 
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literature  as  the  direct  duty  cycle  MPPT  control.  In  this  scheme, 
the  PI  block  in  Fig.  1  is  eliminated  and  duty  cycle  is  computed 
directly  in  the  MPPT  algorithm.  This  scheme  offers  number  of 
advantages:  (1)  it  simplifies  the  tracking  structure,  (2)  it  reduces 
the  computation  time  and  (3)  no  tuning  effort  is  needed  for  the  PI 
gains.  In  short,  it  replaces  the  sophisticated  MPPT  control  with  a 
more  simplified  structure  while  maintaining  similar  optimal 
results. 


3.  Partial  shading  condition 

Fig.  2  shows  a  PV  array  in  a  typical  series-parallel  configura¬ 
tion.  In  this  example,  the  modules  are  connected  in  strings,  with 
three  modules  per  string.  When  one  of  the  modules  in  the  string 
experiences  less  illumination  due  to  shading,  its  voltage  drops; 
thus  it  behaves  as  a  load  instead  of  a  generator.  A  hot  spot  ensued 
and  typically  a  bypass  diode  is  connected  in  parallel  with  each  PV 
module  to  protect  the  shaded  module  from  being  damaged. 
Additionally,  a  blocking  diode  is  connected  at  the  end  of  each 
string  (combination  of  series  modules  in  one  current  path)  to 
provide  the  protection  against  reverse  current  caused  by  the 
voltage  mismatch  between  the  parallel-connected  strings. 


In  normal  condition,  i.e.,  when  the  solar  insolation  on  the 
entire  PV  array  is  uniform,  as  shown  in  Fig.  2(a),  the  P-V  curve 
exhibits  the  typical  unique  MPP  (curve  1  of  Fig.  2(c)).  During 
partial  shading,  as  the  third  PV  module  being  less  illuminated 
(shown  by  shaded  block  in  Fig.  2(b)),  the  difference  in  insolation 
between  two  modules  activates  the  bypass  diode  of  module  3. 
As  a  result,  two  stairs  current  waveform  is  created  on  the  I-V 
curve.  Consequently,  the  corresponding  P-V  curve  is  character¬ 
ized  by  several  local  peaks  and  one  global  peak  (GP),  as  depicted 
by  curve  2.  Furthermore,  if  the  bypass  diode  is  removed,  the  PV 
array  exhibits  only  a  single  peak  (curve  3);  but  this  is  achieved  at 
the  expense  of  a  significant  reduction  in  power.  Therefore,  in 
general,  the  bypass  diode  is  always  installed  to  improve  the 
power  throughput  of  the  PV  array,  despite  the  complication  that 
arises  during  partial  shading  [6]. 


4.  Categorization  of  MPPT  techniques 

In  this  paper,  the  selected  MPPT  are  systematically  categorized 
into  two  main  groups.  First,  which  is  by  far  the  most  popular,  is 
referred  to  as  the  conventional  MPPT.  Three  main  methods 
(together  with  their  important  variations),  namely  the  Perturb 
and  Observe  (P&O),  Incremental  Conductance  (IC)  and  the  Hill 
Climbing  (HC),  will  be  discussed  in  detail.  The  other  group  is 
based  on  soft  computing  techniques.  They  are  becoming  more 
important  lately  due  to  the  availability  of  vast  and  cheap 
computing  power.  The  techniques  that  fall  into  this  category 
include  the  Fuzzy  Logic  Controller  (FLC),  Artificial  Neural  Network 
(ANN)  and  Evolutionary  Algorithms  (EA).  The  last,  in  particular,  is 
of  interests  due  to  its  natural  suitability  to  adapt  for  multi-modal 
problem  imposed  by  the  partial  shading  condition.  They  shall  be 
discussed  comprehensively  in  this  paper. 

Besides  these,  there  are  other  MPPT  described  in  literature: 
Fractional  Short  Circuit  Current  (FSC)  [7-10],  Fractional  Open 
Circuit  Voltage  (FOV)  [11-14]  and  Ripple  Correlation  Control 
(RCC)  [15-19].  These  MPPTs  have  limited  accuracy  but  they  do 
have  their  own  distinct  advantages,  namely  fewer  sensors  and 
simpler  algorithm;  thus  they  offer  a  reliable  and  lower  cost 
solution  for  certain  applications.  Other  techniques  include  current 
sweep  method  [20],  DC-link  capacitor  drop  control  [21,22],  load 
current  and  load  voltage  minimization  [23-27],  dP/dV  or  dP/dl 


a 


b 
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Fig.  2.  Operation  of  PV  array  (a)  under  uniform  insolation  (b)  under  partial  shading  (c)  the  resulting  I-V  and  P-V  curve  for  (a)  and  (b). 
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feedback  control  [28-32],  linear  current  control  [33],  state  based 
MPPT  [34],  best  fixed  voltage  algorithm  [35],  linear  reoriented 
coordinate  method  [36]  and  slide  control  method  [37],  They  are 
not  as  popular  and  hence  have  been  omitted  for  brevity. 


5.  Conventional  MPPT  under  uniform  insolation 

5.3.  Perturb  and  observe  (P  &  0) 

5.1.1.  Principle  of  operation 

The  working  principle  of  P&O  is  depicted  by  the  flowchart  in 
Fig.  3.  The  algorithm  introduces  a  perturbation  (<P)  in  the  operat¬ 
ing  voltage  and  current  of  PV  array  and  subsequently  the  change 
in  the  operating  power  is  observed.  The  increment  in  operating 
power  implies  that  the  converter  is  getting  closer  to  the  MPP. 
Accordingly,  in  the  next  sampling  cycle,  the  direction  (slope)  of 
perturbation  is  retained  and  the  reference  voltage/current  is 
further  increased  by  an  amount  <P.  Note  that  once  the  vicinity  of 
MPP  is  reached,  with  each  new  perturbation  (which  is  changed  in 
sign  alternately),  the  algorithm  will  go  back  and  forth  around  the 
MPP.  Consequently,  it  will  never  stick  exactly  at  MPP;  rather,  it 
oscillates  around  that  point  indefinitely. 


using  a  simple  auxiliary  resonant  circuit,  an  MPPT  for  soft- 
switching  boost  converter  is  proposed.  To  track  the  MPP  under 
rapid  changes  in  insolation,  authors  in  [42]  propose  a  three-point 
weight  comparison  P&O  method.  The  slope  of  the  <P>  is  decided 
based  on  the  comparison  of  the  actual  operating  power  to  two 
preceding  ones.  Other  efforts  to  solve  this  issue  are  carried  out  by 
exploiting  the  sampling  rate  of  MPPT.  For  instance,  the  sampling 
rate  is  optimized  in  [43  j,  while  authors  in  [44],  simply  use  a  high 
sampling  rate. 

Work  in  [45]  suggests  an  adaptive  P&O  scheme  for  a  grid- 
connected  three-phase  inverter.  In  the  beginning,  the  is  set  to 
be  10%  of  the  V0o  Each  successive  perturbation  is  set  to  50%  of  the 
preceding  one  until  the  value  of  <P  is  0.5%  of  the  V0c-  Although  the 
method  exhibits  somewhat  better  performance,  it  is  still  not  fully 
adaptive  due  to  the  predetermined  perturbation  steps.  Patel  and 
Agarwal  [46]  proposed  a  threshold-based  MPPT  scheme  for 
single-stage  grid-connected  PV  system.  The  system  operates  in 
the  continuous  conduction  mode  and  employs  only  voltage 
sensor.  The  array  current  is  estimated  using  the  inductor  current. 


5.2.  Incremental  conductance  (IC) 


5.1.2.  Previous  important  work  on  P&O 

Since  <P  can  be  constant  or  variable,  two  types  of  P&O  are 
reported  in  the  literature,  i.e.,  the  fixed  step  and  adaptive, 
respectively.  For  the  fixed  step  P&O,  authors  in  [38]  proposed  a 
two-stage  grid  connected  inverter;  the  first  stage  is  comprised  of 
a  two-switch  buck-boost  converter  that  performs  the  MPPT  and 
sinusoidal  waveform  generation.  Hashimoto  et  al.  [39]  proposed  a 
novel  MPPT  control  algorithm  for  a  half  bridge  inverter  with  a 
control  unit  consisting  of  a  multistage  buck-boost  converter.  The 
proposed  MPPT  enables  each  of  the  PV  modules  to  generate  its 
maximum  power  by  simply  detecting  only  the  total  output  power 
of  the  PV  system.  Authors  in  [40]  exploit  the  capability  of  multi¬ 
objective  optimization  technique  to  design  the  one-cycle  con¬ 
troller  for  single-stage  inverter.  The  optimization  technique  is  the 
basis  to  design  the  P&O  based  method.  In  another  work  [41], 


Fig.  3.  Flow  chart  of  conventional  Perturb  and  Observe  (P&O)  method. 


5.2.3.  Principle  of  operation 

The  idea  behind  1C  is  to  incrementally  compare  the  ratio  of 
derivative  of  conductance  with  the  instantaneous  conductance 
[47].  It  is  derived  from  the  fact  that  at  MPP,  the  derivative  of 
power  with  respect  to  voltage  (dP/dV)  is  in  fact  zero,  i.e., 


dP 

dV 


d(Vl) 

dV 


-l+v§-° 


0) 


Eq.  (1)  can  be  rearranged  in  the  following  form: 

I  _  dl  ~  A3 
V  ~  dV  =  AV 


(2) 


where  A3  and  AV  are  the  increments  of  PV  current  and  voltage, 
respectively.  The  basic  rules  for  1C  can  thus  be  derived  from  the 
P-V  characteristics  shown  in  Fig.  4  and  can  be  written  as: 

Left  of  MPP  (3) 

[  ^  ^ ,  Right  of  MPP 


Using  the  rules  in  (3),  the  basic  flow  chart  for  1C  is  depicted  in 
Fig.  5.  It  can  be  noticed  that  the  reference  signal  is  based  on 
voltage,  V*.  Since  the  rules  in  (3)  are  derived  using  P-V  curve,  the 
current  cannot  be  used  as  the  final  output.  Instead,  the  P-3 
characteristics  curve  is  utilised.  To  use  IC  as  a  current  based 
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Fig.  5.  The  basic  flow  chart  of  the  voltage  based  IC  method. 


MPPT,  Eq.  (1)  can  be  modified  as: 

^  =  ^>=V+/^=0 
dl  dl  dl 

Accordingly,  the  rules  in  (3)  will  be  transformed  to 
(W=-\,  At  MPP 

Left  of  MPP 
Right  of  MPP 


(4) 


(5) 


Using  (5),  the  flow  chart  in  Fig.  5  can  be  transformed  into 
current  based  IC  by  interchanging  V  with  /,  AV  with  A /  and  V(k)' 
with  /(/<)*.  Subsequently,  the  method  will  be  characterized  by 
incremental  resistance  (1R),  not  IC. 


5.2.2.  Previous  work  on  IC 

Considering  its  superior  dynamic  performance,  many  authors 
[48-53]  employ  IC  in  its  original  form.  For  instance,  in  [48]  an 
interface  circuit  using  IC  for  solar  cells  without  battery  storage  is 
proposed.  In  two  separate  works  [49,50],  the  IC  method  is  applied 
to  track  the  MPP  in  a  single-stage  grid  inverter  system  using  a 
two-mode  scheme.  If  the  insolation  is  sufficient,  it  operates  in 
solar  generation  mode,  otherwise  in  active  power  line  condition¬ 
ing  mode.  In  [51  ],  a  hybrid  approach  is  proposed  to  track  the  MPP 
of  the  module-integrated  converter.  The  authors  utilize  a  linear 
function  to  segmentize  the  I-V  curves  into  MPPs  and  non-MPPs 
region. 

Another  effective  way  to  utilize  the  IC  is  to  generate  an  error, 
e,  using  the  instantaneous  conductance  and  the  incremental 
conductance  [52].  Mathematically,  it  can  be  written  as: 


dl  l 

dV+  V 


=  e 


(6) 


From  (6),  it  can  be  seen  that  the  value  of  e  is  zero  at  MPP. 
The  value  of  e  is  usually  selected  based  on  the  trade-off  between 
the  accuracy  of  (3)  and  the  allowable  oscillations  around  MPP. 
Considering  the  advantages  of  (6),  authors  in  [52]  introduce 


marginal  error  e  in  (3),  i.e„  (d//dV+l/V)  <  e.  Furthermore,  in  [53], 
a  totally  different  approach  is  proposed;  1C  is  implemented  using 
the  mixed-signal  circuit.  A  major  portion  of  the  MPPT  algorithm  is 
realized  using  analog  circuits  which  exploit  the  linear  region  of  a 
field-effect  transistor. 

For  adaptive  IC,  numerous  works  have  been  reported  [54-57]. 
For  instance,  Mei  et  al.  [54]  suggest  an  improved  IR  method, 
where  two  modes  of  perturbation  step  are  used,  i.e.,  fixed  and 
variable.  In  another  work  related  to  adaptive  IC,  Li  and  Wang  [55], 
vary  the  value  of  <f>  in  original  IC  method  with  respect  to  dP/dV 
using  the  following  equality 


&  =  N 


P{k)-P(k- 1) 
V(k)-V{k~I) 


(7) 


In  (7),  N  is  the  scaling  factor,  which  is  tuned  using  a  trial  and 
error  method.  Authors  in  [56]  propose  a  three-phase  grid-con¬ 
nected  PV  system  with  modified  IC  method  and  reactive  power 
compensation.  Similar  to  [55],  a  threshold  value,  based  on  the 
change  in  solar  insolation,  is  used  to  adjust  <t>.  In  the  case  of 
increasing  insolation  i.e.,  AP  is  positive,  <P  is  selected  to  be  smaller 
and  consequently  a  smooth  tracking  towards  MPP  is  achieved.  By 
investigating  the  characteristics  of  PV  array  in  the  voltage  and 
current  source  region,  authors  in  [57]  propose  another  type  of 
adaptive  IC.  If  r=V(k)//(k)  and  rs=dV{k)/dl{k),  then  (8)  outlines  the 
basic  idea  of  their  method,  i.e., 


!-oo  <  r/rs  <  -1,  Voltagesourceregion 
r/rs  =  — 1,  At  MPP  (8) 

-1  <  r/rs  <  0,  Currentsourceregion 


Using  (8),  the  reference  value  of  the  MPPT  algorithm  is  expres¬ 
sed  by: 

rre/(fc+l)  =  rref(k)+Mur  (9) 


M  =  1  +r/rs  (10) 

where  rref  is  a  reference  value  of  a  constant  resistance  and  ur  is  a 
unity  conversion  factor  for  the  resistance  domain. 


5.3.  Hill  climbing  (HC) 


5.3.1.  Operation  of  HC 

In  principle,  the  HC  operation  is  very  similar  to  PSiO;  the  only 
difference  is,  instead  of  perturbing  the  voltage  or  current,  it 
updates  the  operating  point  of  the  PV  array  by  perturbing  the 
duty  cycle,  d(k )  by  a  fixed  step-size  ($)  in  the  direction  of 
increasing  power.  The  perturbation  direction  is  reversed  if 
P(k)  <  P(k—  I);  an  indication  that  the  tracking  is  not  moving 
towards  the  MPP.  This  can  be  described  by  the  following  equa¬ 
tion: 

j  d(k)+<P.  if  P(k)  >  P(k-i) 

°(k)~  j  d(k)-<P.  if  P(k)  <  P(k-I)  (11) 

Due  to  the  fact  that  the  duty  cycle  is  changed  directly  without 
a  PI  controller,  the  HC  method  is  sometimes  known  as  the  direct 
duty  cycle  technique.  As  with  the  case  of  PSiO,  this  method  can  be 
characterized  by  fixed  or  adaptive  perturbation  step. 


5.3.2.  Previous  work  on  HC 

Considering  the  fact  that  the  parallel  capacitance  connected  at 
the  output  affects  the  characteristics  of  the  PV  array,  Kasa  et  al. 
[58]  estimate  <P>  by  the  following  equation: 

<P  =  Km  x  Cpy  x  - -  (12) 

Cpv0<) 

where  CPV(k)  is  the  estimated  value  of  capacitance  CPV.  It  is 
reported  that  the  power  from  the  PV  array  can  be  transferred  to 
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the  utility  grid  with  the  nearly  unity  power  factor  at  the 
maximum  power  point.  Koutroulis  et  al.  [59]  use  HC  for  buck 
converter  battery  charging  PV  system,  while  in  [SO],  the  authors 
simplify  it  for  a  fly  back  inverter  using  only  voltage  sensor.  Jain 
and  Agarwal  [61  ]  apply  the  conventional  HC  scheme  on  a  high 
gain  inverter,  operating  in  discontinuous  current  mode.  This 
allows  for  the  low  PV  array  voltages  (typically  50-100  V  range) 
to  be  boosted  up  to  the  level  of  grid  voltage  (115  V  or  230  V  ac 
systems).  In  another  work  [62],  HC  is  used  in  a  parallel  connected 
MPPT  system  for  stand-alone  PV  power  generation.  The  proposed 
structure  reduces  the  negative  influence  of  dc-dc  converter 
losses,  hence  increasing  the  efficiency  of  power  generation.  Kwon 
et  al.  [63]  propose  a  simple  three-level  boosting  MPPT  control  for 
a  three  phase  grid  connected  PV  system.  The  proposed  MPPT 
converter,  which  utilizes  the  HC  algorithm,  reduces  the  reverse 
recovery  losses  of  the  diodes  and  thereby  increases  the  overall 
efficiency. 

Chiang  et  al.  [64]  propose  a  modification  to  (12)  to  formulate 
an  adaptive  HC  algorithm.  The  perturbation  step  is  computed 
using  the  following  relationship: 


<£  = 


AP/Ad 

P/d 


(13) 


In  a  separate  work,  Xiao  and  Dunford  [65]  exploit  the  P-D 
curve  for  the  PV  array.  Using  the  value  of  (dP/dD),  they  tuned  the 
perturbation  value  <F>  using 

0<k>  =  Mwl^)  (14> 


If  the  change  in  array  power  (AP)  is  large  due  to  rapid 
environmental  variations,  <P(k)  is  forced  likewise  to  ensure 
improved  dynamic  performance.  When  AP  is  small,  the  system 
is  near  to  MPP  and  accordingly,  the  value  of  <P(k)  will  also  be 
small.  Wolfs  and  Tang  [66]  propose  another  modification  to  (14) 
using  the  relation  of  dP/dV  as  follows: 


<P  =  M 


AP 

AV 


(15) 


However,  since  the  value  of  M  in  (15)  primarily  decides  the 
performance  of  MPPT,  manual  tuning  of  this  scaling  factor  can  be 
extremely  sensitive  to  initial  operating  conditions.  To  overcome 
this  issues,  Pandey  et  al.  [67]  propose  a  fully  adaptive  HC  method 
to  automatically  tune  M  under  all  operating  conditions.  This  is 
calculated  using 

a /t _  \avmax\&Dmax  1fi 

|APMax|  (  ’ 


where,  AVMax  and  APMax  are  the  maximum  change  in  array 
voltage  and  power,  respectively,  with  respect  to  maximum 
change  in  duty  cycle  ADmax- 

As  the  HC  method  works  on  the  same  principle  like  P&O,  its 
dynamic  performance  suffers  similar  problem  to  the  latter,  i.e., 
the  divergence  of  MPP.  The  three  point  sampling  technique  [42] 
can  also  be  applied  in  HC  to  overcome  this  problem. 


5.4.  IC  with  direct  duty  cycle 

Another  class  of  direct  duty  cycle  control  is  inspired  by  the  IC 
method.  The  duty  cycle  of  the  converter  is  varied  according  to  the 
rules  given  in  (3).  Like  other  techniques,  this  method  can  also  be 
employed  in  its  original  form  or  adaptive.  Xuesong  et  al.  [68] 
utilize  the  IC  algorithm  depicted  by  the  flowchart  in  Fig.  5. 
Despite  the  absence  of  PI  control  loop,  it  is  shown  through 
simulations  that  this  method  can  satisfactorily  track  the  MPP. 
Using  the  same  concept  of  marginal  error  e  as  reported  in  [52], 
Safari  and  Saad  [69]  use  a  Cuk  converter  as  the  MPPT  converter. 


Although  fast  and  accurate,  the  method  is  not  generic  because  for 
different  PV  system,  the  perturbation  step  <t>  and  marginal  error  e 
need  to  be  recomputed  based  on  trial  and  error  approach. 

Various  adaptive  IC  are  proposed  to  obtain  the  best  trade-off 
between  MPP  oscillations  and  dynamic  performance.  In  [70], 
based  on  a  threshold  value  of  operating  power  e,  a  two  step 
based  direct  duty  cycle  1C  method  is  proposed.  Small  value  of  <P  is 
used  in  the  vicinity  of  the  MPP  while  big  step  is  used  in  other 
regions.  Using  the  same  idea  of  [67],  Liu  et  al.  [71]  propose 
another  auto  tuning  scheme  for  the  scaling  factor  M  in  (17). 
The  value  of  M  is  constrained  to  the  following  inequality: 


M  <  A Dmax/ 


dP 

dV 


Fixed  Step  =  ADmax 


(17) 


If  (17)  is  not  satisfied,  the  algorithm  works  with  a  fixed  step 
size.  Additionally,  for  smooth  start,  a  simple  start  program  is 
introduced  to  the  MPPT  method.  Another  work  carried  out  by  [72] 
presents  an  alternative  form  of  adaptive  IC;  <P  is  calculated  using 


<P=  +M 


Al  I(k ) 
A V+V(k) 


=  ±Cf(k ) 


(18) 


As  duty  cycle  always  fall  in  a  particular  range,  M  is  allowed  to 
be  selected  from  the  following  inequalities: 


ADmin  m  A  DMax 
/max(^)  f  MAx(fc) 

where 


(19) 


/  MAX 


A  /  1(1) 
AV  V(l) 


(20) 


6.  Soft  computing  MPPT  under  uniform  insolation 

6.1.  Fuzzy  logic  control  (FLC) 


6.1 .1.  Operation 

In  FLC,  a  mathematical  model  of  the  system  is  not  required. 
This  in  itself  is  a  significant  advantage  because  the  uncertainties 
such  as  un-modelled  physical  quantities,  non-linearity  and  unpre¬ 
dictable  changes  in  operating  point  can  be  excellently  dealt  with 
[73,74].  However,  it  does  require  the  designer  to  have  some  prior 
knowledge  of  how  the  output  responds  qualitatively  to  the 
inputs.  The  typical  process  structure  of  an  FLC  is  shown  in 
Fig.  6.  The  controller  has  three  functional  blocks  namely  fuzzifica¬ 
tion,  rules  inferences  and  defuzzification.  In  addition,  it  has  a  rule 
table  in  which  the  designed  rules  are  stored.  The  process  in  which 
the  FLC  performs  the  calculation  is  called  rule  inference. 

The  inputs  to  a  FLC-based  MPPT  are  usually  an  error  E  and  a 
change  in  error  A£.  Since  dP/dV  vanishes  at  the  MPP,  both  inputs 
can  be  calculated  as  follows: 


E(k) 


P(k)-P(k-t) 
V(fc)-V(k- 1) 


or  E(k) 


P(k)-P(k- 1) 
Kk)-I(k-l) 


(21) 


AE  =  E(k)—E(k—i)  (22) 

Once  E  and  A£  are  computed  and  converted  to  the  linguistic 
variables,  the  FLC  output  can  also  be  defined  linguistically  in 
terms  of  voltage  (AV),  current  (Al)  or  duty  cycle  (AD).  The 


Fig.  6.  Fuzzy  logic  control  (FLC)  process  structure. 
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linguistic  variables  for  the  output  signal  are  usually  assigned 
based  on  innate  knowledge  of  the  particular  PV  system  being 
employed. 

6.1.2.  Previous  work  using  FLC 

in  [75],  an  FLC  with  seven  linguistic  variables  is  used,  while  in 
[76]  and  [77],  the  membership  functions  are  made  less  symmetric 
to  give  more  significance  to  the  specific  linguistic  variables.  In  an 
effort  to  overcome  tracking  deterioration  with  respect  to  changes 
in  solar  insolation  experienced  by  [78],  authors  in  [79]  introduce 
the  array  power  variation  and  duty  cycle  as  the  inputs  for  the  FLC. 
Although,  it  performs  satisfactorily  in  varying  environmental 
conditions,  it  could  not  reduce  the  steady  state  error  in  the  array 
power.  Extending  the  idea  of  [78]  and  [79],  authors  in  [80] 
propose  a  3-input  FLC  based  MPPT.  The  information  derived  from 
array  current,  power  and  duty  cycle  of  converter  are  selected  as 
the  inputs  to  the  FLC.  Some  improvements  are  seen  in  [81  ],  where 
the  authors  employ  the  Fuzzy  Cognitive  Networks  (FCN)  to 
enhance  the  performance  of  [78].  The  voltage,  current,  tempera¬ 
ture,  solar  insolation  and  control  variable  (current)  are  used  as  the 
nodes  of  the  FCN.  The  weights  of  the  interconnection  nodes  are 
determined  using  a  wide  range  of  different  climatic  conditions. 
Once  the  FCN  is  trained,  it  can  be  mounted  on  any  PV  system. 
The  proposed  method  results  in  better  tracking  speed,  but  at  the 
expense  of  additional  switch  and  a  sensor. 

In  [82],  the  authors  exploit  the  advantages  of  the  FLC  to  improve 
the  performance  of  the  conventional  HC  method.  The  FLC  is  devel¬ 
oped  by  translating  the  HC  algorithm  into  16  fuzzy  rules.  In  addition, 
several  authors  utilise  the  idea  of  adaptive  FLC.  For  instance,  in  [83], 
an  adaptive  FLC  is  proposed  to  facilitate  the  constant  tuning  of  the 
membership  functions  and  the  rule  based  table  in  order  to  achieve 
optimum  performance.  In  another  work,  [84],  the  scaling  factor  of 
both  fuzzy  inputs  and  output  are  automatically  tuned  to  achieve  the 
better  dynamic  performance  of  MPPT. 

Authors  in  [85]  propose  a  single  input  fuzzy  logic  controller 
(SI-FLC)  by  applying  the  “signed  distance  method"  [86].  The  input 
to  SI-FLC  consists  of  only  one  variable  known  as  "distance”.  This  is 
in  contrast  to  the  conventional  FLC,  which  requires  an  error  and 
the  derivative  (change)  of  the  error  as  its  inputs.  The  reduction  in 
the  number  of  inputs  simplifies  the  rule  table  to  one-dimensional, 
allowing  it  to  be  treated  as  a  single-input  single-output  controller. 

6.2.  Artificial  neural  network  (ANN) 

6.2.1.  ANN  operation 

The  ability  of  ANN  to  recognize  and  estimate  unknown  para¬ 
meters  inspired  its  application  for  MPP  tracking  [87].  The  input 
variables  to  the  ANN  can  be  PV  array  parameters  such  as  V0c  and 
/sc.  insolation  and  temperature,  or  any  combination  of  these. 
The  output  is  usually  a  reference  signal,  which  can  be  either  the 
voltage,  current  or  duty  cycle.  To  accurately  identify  the  MPP,  the 
weights  associated  to  the  neurons  have  to  be  carefully  computed 
through  a  comprehensive  training  process.  Once  trained,  the  ANN 
can  be  used  as  the  MPP  estimator  which  will  give  the  reference 
value  ( VMP  or  IMP )  to  the  MPPT  controller. 

6.2.2.  Previous  work 

Typically,  the  ANN  is  used  to  estimate  the  MPP  with  respect  to 
environmental  variations.  For  example,  in  [88],  the  ANN  is 
utilized  to  identify  the  MPP  using  a  gradient  descent  algorithm 
training.  Other  authors  used  the  capability  of  ANN  to  improve  the 
P&O  or  IC  method;  Alabedin  et  al.  [89]  develop  an  ANN  based  P&O 
controller.  Compared  to  the  FLC  based  MPPT,  the  proposed 
method  results  in  improved  performance  in  dealing  with  the 
fluctuations  in  the  array  power.  Authors  in  [90]  propose  a  novel 


Layer  h  Layer  j  Layer  k 


ANN  technique,  as  shown  in  Fig.  9,  to  further  improve  the  P&O 
performance.  It  can  be  seen  that  the  insolation  and  temperature 
are  used  as  the  input  for  the  network  together  with  20  neurons  in 
the  hidden  layer.  In  another  work,  [91  ],  the  authors  improve  the 
conventional  IC  method  by  combining  it  with  a  trained  ANN. 

Most  ANN  methods  [88-91]  employ  back-propagation  (BP) 
training  algorithm.  Larger  number  of  hidden  nodes  yield  more 
accurate  results,  but  at  the  expense  of  longer  computational  time. 
Hence,  in  rapid  fluctuations  of  environment,  the  ANN  MPPT  may  not 
be  able  to  respond  quickly  enough  to  cope  with  the  fast  changes. 
Moreover,  the  BP  algorithm  remembers  the  new  value  but  forgets 
the  old,  leading  to  the  poor  memory  retainement.  Therefore,  in 
general,  ANN  is  unsuitable  for  low  cost  microprocessors  Fig.  7. 

6.3.  Evolutionary  algorithm  (EA ) 

Another  soft  computing  approach  that  is  gaining  attention 
recently  is  the  evolutionary  algorithm  (EA).  It  is  a  stochastic 
method  that  appears  to  be  very  efficient  in  optimizing  real-valued 
non-linear  and  multi-modal  objective  functions  [92,93).  Since  the 
technique  is  based  on  search  optimization,  in  principle,  it  should 
be  able  to  locate  the  MPP  regardless  of  environmental  variations. 
Various  EA  methods  for  MPPT  are  found  in  the  literature;  the 
most  popular  ones  are  Particle  Swarm  Optimization  (PSO)  [94], 
Genetic  Algorithm  (GA)  and  Differential  Evolution  (DE) 

6.3.1.  PSO 

PSO  is  a  population-based  search  method,  modelled  after  the 
behaviour  of  bird  flocks  [95],  The  algorithm  maintains  a  swarm  of 
individuals  (called  particles),  where  each  particle  represents  a 
candidate  solution.  Particles  follow  a  simple  behaviour:  emulate 
the  success  of  neighbouring  particles,  and  achieve  their  own 
successes.  The  position  of  a  particle  is  therefore  influenced  by 
the  best  particle  in  a  neighbourhood,  as  well  as  the  best  solution 
found  by  the  particle.  Particles  positions,  xb  are  adjusted  using: 

x,-k+1  =  xk+vk+1  (23) 

Note  that  in  (23)  where  the  velocity  component,  v b  represents 
the  step  size.  The  velocity  is  calculated  by: 

vk+ 1  =  wvk  +  Ci  ri  [P(,est  j-xf)  +  c2r2{Gbest-xk )  (24) 

where  w  is  the  inertia  weight,  c ^  and  c2  are  the  acceleration 
coefficients,  r,.  r2eU(0,  1),  Pbesti  is  the  personal  best  position  of 
particle  i,  and  Gbest  is  the  neighbourhood  best  position  of  particle  i. 
The  inertia  weight  w  plays  an  important  role  in  balancing  the  global 
search  and  local  search.  A  large  w  facilitates  a  global  search  while  a 
small  inertia  weight  improves  a  local  search.  It  can  be  a  positive 
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Fig.  8.  The  basic  flow  chart  of  PSO  method. 


Fig.  9.  The  1-V  and  P-V  curves  used  for  evaluating  the  performance  of  various 
MPPT  methods. 


constant  or  a  positive  decreasing  linear  function  of  iteration  index  j. 
Fig.  8  depicts  the  basic  flow  chart  of  the  PSO  method. 

Authors  in  [96]  employ  the  PSO  algorithm  to  predict  the  MPP 
of  PV  system.  At  the  beginning  of  the  evolution,  the  conventional 
PSO  is  applied  for  the  global  search.  Once  it  arrives  at  the  vicinity 
of  the  MPP,  a  special  algorithm  is  triggered  to  quickly  obtain  the 
local  optimal  point.  Such  scheme  effectively  improves  the  speed 
of  local  search.  In  general,  PSO  is  more  suitably  used  for  partial 
shading  conditions,  as  shall  be  described  in  the  following  sections. 


6.3.2.  Other  EA  techniques 

In  [97],  a  hybrid  CA-ANN  MPPT  is  proposed.  The  GA  is  used  to 
obtain  the  optimized  values  of  array  power  and  voltage  for 
different  insolation  and  temperatures  conditions.  The  results  are 
used  for  the  offline  training  of  the  ANN.  In  a  similar  manner, 
authors  in  [98]  propose  another  GA-ANN  based  MPPT,  while  in 
[99],  GA  is  used  to  optimize  the  FLC  based  MPPT. 


7.  MPPT  during  partial  shading 

In  this  section,  the  performance  conventional  MPPT  (P8iO,  IC, 
and  HC)  and  the  soft  computing  techniques  (FLC,  ANN  and  EA) 
under  partial  shading  are  reviewed. 


7.1.  Modified  P&O 


To  aid  the  discussion,  Fig.  9  is  used.  In  this  figure,  the  1-V  and 
P-V  curves  are  shown:  (1)  for  uniform  insolation  (solid  trace),  and 
(2)  for  partially  shaded  condition  (dotted  trace).  Under  uniform 
insolation  only  one  peak  exists,  i.e.,  GP1 ;  the  conventional  P&O 
almost  always  successfully  track  this  point.  For  partial  shading, 
five  peaks  labelled  by  A-E  emerge;  clearly,  peak  E  is  the  true  MPP 
(GP)  while  the  others  are  local.  When  shading  occurs,  the  array 
operating  point  shifts  from  1  to  2;  the  perturbation  sign  for  the 
P&O  reverses  and  control  is  increased.  Eventually  it  will  recognize 
point  D  as  the  GP,  which  is  obviously  wrong. 

The  above  discussion  highlights  the  limitation  of  the  conven¬ 
tional  P&O  during  partial  shading  and  hence  the  need  for  remedial 
effort.  In  [100]  a  two-mode  modified  P&O,  outlined  by  the  flowchart 
in  Fig.  10  is  proposed.  The  complete  algorithm  is  divided  into  two 
parts:  (1)  the  main  program  and  (2)  the  GP  tracking  routine. 
To  avoid  scanning  the  whole  P-V  curve,  two  boundary  conditions 
(Vmln  and  Vmax)  are  introduced;  these  values  are  based  on  the 
voltage  in  the  P-V  curve  of  the  deployed  PV  array.  Despite  this 
simplification,  for  certain  shading  conditions,  the  algorithm  still 
needs  to  scan  almost  80%  of  the  P-V  curve  to  ensure  none  of  the 
potential  peak  GP  is  missed.  As  a  result,  the  tracking  speed  is  greatly 
compromised.  The  scanning  of  P-V  cun/e  is  terminated  when  the 
highest  local  peak,  i.e.,  the  true  GP,  is  found. 

In  another  work  [101],  an  alternative  P&O  using  the  compar¬ 
ison  of  two  instantaneous  power  values  is  proposed.  Mathemati¬ 
cally,  the  comparison  can  be  written  as 


Pm(t)-Pref(t) 
Pm(t~  1) 


(25) 


In  (25),  [Pm  (t)[  and  [Pref  (t)[  are  the  instantaneous  measured 
power  and  instantaneous  maximum  power  reference,  respec¬ 
tively.  Despite  the  good  GP  tracking  capability,  it  introduces 
various  new  coefficients  (which  are  not  shown  here  for  simpli¬ 
city)  that  complicate  the  overall  MPPT  process.  Moreover,  the 
method  is  tested  only  on  one  shading  condition;  thus,  its  effec¬ 
tiveness  is  not  fully  validated.  Authors  in  [102]  propose  a  voltage 
sweep  method  to  periodically  change  the  voltage  of  PV  array  from 
its  maximum  value  (near  to  V0c)  to  a  minimum  value  (near  to  /sc). 
At  each  sampling  cycle,  the  operating  voltage  and  current  are 
measured  and  stored  in  the  microcontroller  memory.  After 
identifying  the  region  of  global  maxima,  the  P&O  method  is 
executed  to  maintain  the  operation  at  the  GP. 


7.2.  Modified  IC 

Since  the  global  and  local  peaks  exhibit  the  same  derivative 
characteristics,  i.e.,  ( dP/dV=0 )  or  (dP/d/= 0),  the  original  IC  algo¬ 
rithm  is  unable  to  recognize  the  true  MPP.  Accordingly,  several 
authors  propose  modified  IC  schemes  during  partial  shading.  For 
example,  Kobayashi  et  al.  [103]  implement  a  two-stage  IC 
method.  In  the  first  stage,  the  PV  array  is  forced  to  operate  into 
the  neighbourhood  of  the  GP  using  the  values  of  the  maximum 
voltage  and  current,  i.e., 

RMp  =  k ^  (26) 

Imp 

where  k  is  the  correction  factor,  while  VMP  and  IMP  are  approxi¬ 
mated  to  be  80%  of  V0c  and  90%  of  ISc,  respectively.  Subsequently, 
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Main  program 


Fig.  10.  Flow  chart  to  track  the  GP  under  partial  shading  condition  using  modified  P&O  [100], 


in  the  second  stage,  adjustment  is  made  to  move  the  operating 
point  towards  the  GP  based  on  the  IC  method.  In  another  attempt, 
authors  in  [104]  propose  a  linear  function  to  track  the  GP  under 
partial  shading.  Mathematically,  it  can  be  written  as: 

V*  =  V.grid  x  /(/<)  (27) 

*out 


7.3.  Modified  HC 

Since  the  HC  method  is  based  on  the  P&O  principle,  the 
algorithm  also  traps  at  the  local  minima.  In  [105],  Lei  et  al. 
proposed  a  duty  cycle  sweep  method  based  on  modified  HC  to 
cater  for  partial  shading  condition.  The  initial  value  of  the  duty 
cycle  is  computed  using  the  following  equation: 


where  Vgrid  and  Iout  is  the  output  voltage  and  current  of  the  grid, 
respectively.  The  linear  function  is  activated  whenever  shading 
condition  occurs.  It  is  detected  using  the  following  relationships: 

V(k)—V(k—1)  <  Vthr  (28) 


I(k)-I(k-r) 
I(k- 1) 


<h, 


(29) 


Once  the  PV  array  is  subjected  to  partial  shading,  (27)  shifts 
the  operating  point  to  the  reference  value  (V).  Then,  the  IC 
method  is  employed  to  track  the  GP.  Although  quite  effective, 
the  method  can  only  be  applied  to  a  grid-connected  system. 
Interestingly,  it  is  observed  in  [104]  that  under  rapid  fluctuation 
of  insolation,  the  performance  of  the  modified  IC  deteriorates 
more  rapidly  compared  to  the  conventional  method. 


(30) 


where  RLoad  is  estimated  using  the  rating  of  PV  array.  The  sweep 
duty  cycle  is  carried  out  for  the  range  of  0-90%.  Later,  the 
conventional  HC  is  improved  by  dividing  the  duty  cycle  perturba¬ 
tion  (cP)  by  two  whenever  the  direction  of  the  <1>  is  increased. 
Similar  to  [100],  this  method  also  needs  to  scan  over  80%  of  the  PV 
curve.  In  another  work,  authors  in  [106]  implemented  a  multiple- 
input  boost  converter  for  micro-inverters  based  on  modified  HC. 


7.4.  ANN  and  FLC  in  partial  shading 

By  carefully  analyzing  (26),  it  can  be  deduced  that  the  operat¬ 
ing  principle  of  FLC  is  identical  to  the  IC;  they  both  work  on  the 
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basis  of  dP/dV  information.  A  zero  value  of  £(/<)  implies  that  the 
algorithm  is  at  the  MPP,  i.e.,  dP/dV=  0.  However,  in  partial 
shading,  all  peaks  (local  and  GP)  exhibit  the  same  dP/dV  char¬ 
acteristic.  Hence,  FLC  too  cannot  guarantee  the  true  GP. 

Accordingly,  a  modified  FLC  based  on  three  stages  operation  is 
proposed  in  [107].  The  key  steps  involved:  (1)  the  scanning  of  P-V 
curve  (in  uniform  or  partial  shading),  (2)  storing  of  all  possible 
local  peaks  and  (3)  performing  the  perturbation  and  observation 
of  the  array  power.  Three  variables,  given  by  the  following 
equations,  are  fed  to  the  FLC  inputs: 

AP  =  P(k)-P(k-l)  (31) 

A/  =  /(/<)— /(fe—1)  (32) 

APm  =  PM(k)-P(k- 1)  (33) 

where  AP  and  A /  are  the  changes  in  array  power  and  current  in 
two  subsequent  samples,  respectively.  A PM  is  the  difference 
between  the  stored  GP  (PM)  and  the  currently  measured  power. 
Using  these  inputs,  thirty-four  FLC  rules  for  the  duty  cycle  signal 
(AD)  are  formulated  based  on  the  HC  algorithm. 

In  addition,  efforts  have  been  made  to  combine  the  conven¬ 
tional  MPPT  with  ANN.  Syafarudin  et  al.  [108]  propose  a  three 
layer  feed-forward  ANN,  combined  with  an  FLC.  The  method  is 
tested  on  different  solar  cell  technologies,  namely  mono¬ 
crystalline  silicon,  thin-film  and  amorphous  silicon.  In  another 
work  by  the  same  authors  [109],  a  Fuzzy  wavelet  network  to 
identify  the  GP  of  non-crystalline  Si  PV  modules  is  proposed.  A 
data  set  which  consists  of  200  input-output  pairs  is  used  to  train 
and  test  the  proposed  ANN  network.  However,  for  both  works,  it 
is  clearly  recognized  that  obtaining  the  reliable  training  data  set 
for  partial  shading  is  the  most  challenging  tasks.  Although 
shading  due  to  towers,  trees  and  roofs  can  be  modelled  easily, 
the  conditions  due  to  the  uncertainty  of  passing  clouds  are 
extremely  difficult  to  estimate.  Since  the  training  is  carried  out 
using  predefined  shading  scenarios,  the  effectiveness  of  the  ANN 
algorithm  could  not  be  guaranteed. 

7.5.  PSO 


Fig.  11.  Bayesian  network  for  MPPT  information  fusion. 


a  Bayesian  fusion  technique  to  combine  PSO  with  1C  methods,  as 
shown  in  Fig.  11.  The  network  has  N=  20  nodes  represented  by 

[Ai  A2 . A,v].  They  are  divided  into  two  sectors,  left  and  right.  Ten 

nodes  [A1?  A2 . A]0]  correspond  to  10  points  of  the  l-V  curve, 

estimated  by  the  PSO  method.  The  remaining  10  nodes  refer  to 
the  10  points,  estimated  by  the  1C  method.  Finally,  the  symbol  S 
denotes  the  estimated  location  after  fusion  process.  Despite  its 
effectiveness,  the  complexity  of  the  method  is  the  major  con¬ 
straint  for  implementing  it  using  a  low-cost  microprocessor. 

Ishaque  et  al.  [118]  adopt  the  PSO  algorithm  in  the  direct 
control  structure,  where  the  duty  cycles  are  used  as  the  positions 
of  the  PSO  particles.  The  main  idea  is  shown  by  flow  chart  in 
Fig.  12.  If  a  drastic  change  in  operating  power  is  observed  (which 
implies  the  occurrence  of  partial  shading),  the  algorithm  begins 
the  search  by  sending  duty  cycles  to  the  power  converter.  The 
information  (voltage  and  current)  obtained  from  these  duty  cycles 
are  used  to  compute  the  Pbest  and  Cbest  values.  Using  the  velocity 
and  position  equations,  the  new  duty  cycles  are  computed  based 
using  the  following  equations: 

$ik+1  =w$ik+Cir1{Pbesti-xik}+c2r2[Gbest-xik}  (36) 


Since  partially  shaded  l-V  curve  exhibits  multimodal  charac¬ 
teristics,  PSO  methods  are  envisaged  very  effective  to  track  the  GP 
under  this  condition.  Miyatake  et  al.  [110,111]  utilize  the  con¬ 
ventional  PSO  to  track  the  GP  in  a  constant  bus  voltage  applica¬ 
tion.  Authors  in  [112]  formulate  an  analytical  expression  of  the 
objective  functions  based  on  PV  current,  insolation  and  tempera¬ 
ture;  then  the  conventional  PSO  is  utilized  to  track  the  MPP.  In 
another  work,  current  based  PSO  method  [113]  is  proposed.  The 
current  through  the  series  inductor  in  boost  converter  is  used  as 
the  reference  signal  to  generate  the  PWM  signals  for  the  switch¬ 
ing  converter.  Phimmasone  et  al.  [114,115]  introduce  a  repulsive 
term  in  the  velocity  equation  of  the  PSO  as  follows: 

vk+ 1  =  wvf  +  Cl  IT  [Pbest  i-Xik}  +  c2r2[Gbest-Xik} 

-c-ir3{centk-Xik}  /  [\centk-Xik  +d}3  (34) 


where  c3  is  the  coefficient  of  repulsive  particle.  The  function 
"cent"  is  the  centre  of  all  particles,  described  by 


centk 


N  vk 


(35) 


Another  work  by  [116]  propose  an  Adaptive  Perceptive  PSO 
(APPSO);  however,  due  to  an  additional  dimensional  search  space, 
the  number  of  particles  is  increased  significantly  compared  to  its 
original  counterpart  [111].  This  increases  the  computation  bur¬ 
den,  thus  making  real-time  implementation  more  difficult. 
Further  improvement  is  seen  in  [117],  where  the  authors  utilize 


dk+l  =  dik+4>ik+1  (37) 

Compared  with  conventional  HC  method,  this  method  is  able 
to  handle  severe  partial  shading  conditions  with  excellent 
dynamic  tracking  speed.  Moreover,  the  steady  state  oscillations 
are  found  to  be  exceptionally  low. 

More  recently,  an  improved  version  of  [118]  (proposed  by  the 
same  authors)  is  published  in  [119].  The  main  idea  is  to  remove 
the  random  number  in  the  accelerations  factor  of  the  conven¬ 
tional  PSO  velocity  equation.  Hence  it  is  appropriately  named  as 
deterministic  the  PSO  (DPSO).  The  modification  is  described  as 
follows: 

Pik+1  =  wv,k  +  [Pbest  i-Xik}  +  [Gbest—xk} 

JJ.  for  0  <v  <V  max  (38) 

Vife+1  =  WVik  +  {Gbest  +  Pbest  i-2Xikj 

The  maximum  change  in  velocity  (Vmax)  is  restricted  to  a 
particular  value,  which  is  determined  based  on  the  critical  study 
of  P-V  characteristics  during  partial  shading.  The  proposed 
method  is  shown  to  have  following  advantages:  (1)  consistent 
solution  is  achieved  even  for  a  small  number  of  particles  (2)  only 
the  inertia  weight  need  to  be  tuned  and  (3)  simpler  MPPT 
structure  compared  to  the  conventional  PSO.  Despite  using  only 
three  particles,  it  is  experimentally  shown  that  the  method 
performs  very  satisfactorily  under  very  severe  shading  conditions. 
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7.6.  Other  EA  methods 

Authors  in  [120]  propose  another  EA  method,  namely  differ¬ 
ential  evolution  (DE).  The  method  is  tested  on  three  P-V  curves, 
each  exhibits  three  peaks.  However,  for  each  curve  presented,  the 
GP  is  on  the  right  side,  which  eventually  turns  out  to  be  only  one 
shading  condition.  Hence,  these  trivial  tests  cannot  be  considered 
sufficient  to  prove  the  effectiveness  of  the  method  if  other  shading 
conditions  occur. 

Using  the  similar  searching  methodology  of  EA  methods, 
another  method  that  exploits  the  capability  of  dividing  rectangles 
(DIRECT)  algorithm  is  developed  in  [121],  It  is  shown  that  the  P-V 
characteristics  curve  can  be  represented  by  the  Lipschitz  func¬ 
tion.  The  algorithm  operates  either  in  the  global  and  local  mode. 
Once  the  occurrence  of  partial  shading  is  confirmed,  the  global 
mode  is  activated  and  DIRECT  algorithm  is  applied  to  locate  the 
vicinity  of  the  GP.  Thereafter,  the  algorithm  switches  to  the  local 
mode  where  a  conventional  P&O  with  a  relatively  small  perturba¬ 
tion  is  employed  to  maintain  the  operating  point  at  the  GP.  Once 
again,  despite  the  improved  tracking  capability,  the  complexity  of 
the  method  is  the  major  issue. 


8.  Summary 

Judging  by  the  amount  of  recent  work,  it  can  be  concluded  that 
the  MPPT  is  still  a  very  active  research  area.  There  are  rooms  for 
improvement,  particularly  when  dealing  with  partial  shading 
conditions.  For  uniform  insolation,  there  is  an  increasing  trend 
in  combining  the  conventional  methods  with  soft  computing 
techniques.  Examples  of  successful  schemes  have  been  demon¬ 
strated  by  various  works  described  in  this  paper.  Such  “hybrid 
MPPT”  is  quite  attractive  due  to  the  availability  of  the  vast 
computational  power  using  relatively  low  cost  microprocessors. 

As  for  the  partial  shading,  it  appears  that  the  EA  methods  are 
the  most  promising.  The  combination  of  direct  duty  cycle  method 
and  EA  can  be  an  interesting  option  to  explore.  The  PSO,  which 
has  a  very  similar  structure  to  HC  has  proven  to  be  a  viable 
option.  However,  there  are  still  potential  areas  of  concerns  for  EA, 


*  A  Lipschitz  function/is  such  that  f(x)  =  f(y) |  <  C|x  =y\  for  all  x  and  y,  where 
C  is  a  constant  independent  of  x  and  y. 


particularly  the  selection  of  appropriate  evolutionary  control 
parameters  such  as  population  size,  mutation  factor  etc.  This 
may  be  resolved  with  the  assistance  of  other  optimization 
technique  such  as  ANN.  There  is  also  a  need  to  consider  the 
viability  of  EA  methods  when  it  comes  to  real  time  implementa¬ 
tion.  The  way  forward  is  to  find  an  effective  method  such  that  the 
search  space  can  be  narrowed,  so  that  the  computational  time  is 
reduced.  Furthermore,  there  are  several  more  EA  techniques, 
which  are  yet  to  be  exploited  fully.  One  example  is  the  differential 
evolution  (DE).  Although  it  has  been  successfully  applied  in  the 
area  of  PV  cell  modelling,  it  has  not  been  used  for  MPPT-except 
for  one  very  basic  work  done  by  [120]. 
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